What does a strongly excited 't Hooft-Polyakov magnetic monopole do? 
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The time evolution of strongly exited SU (2) Bogomolny-Prasad-Sommerfield (BPS) magnetic monopoles in 
Minkowski spacetime is investigated by means of numerical simulations based on the technique of conformal 
compactification and on the use of hyperboloidal initial value problem. It is found that an initially static 
monopole does not radiate the entire energy of the exciting pulse toward future null infinity. Rather, a long- 
lasting quasi-stable 'breathing state' develops in the central region and certain expanding shell structures - 
built up by very high frequency oscillations - are formed in the far away region. 
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The investigation of solitons in particle physics is of fun- 
damental interest (see e.g. [J for a recent review). In par- 
ticular, considerable attention has been paid to the study 
of 't Hooft-Polyakov magnetic monopole solutions of cou- 
pled Yang-Mills-Higgs (YMH) systems 0, Q. However, 
the relevant investigations have been almost exclusively 
restricted to time independent configurations. Thereby it 
is of obvious physical interest to study dynamical proper- 
ties of magnetic monopole configurations. 

This paper is to report about the results of investiga- 
tions concerning the time evolution of a strongly excited 
spherically symmetric SU{2) BPS magnetic monopole 
0, on a fixed Minkowski background spacetime by 
means of numerical techniques. In particular, the under- 
lying Yang-Mills-Higgs system is chosen so that the Yang- 
Mills field is massive while the Higgs field is massless. The 
dynamics starts by hitting the static monopole by a con- 
centrated pulse. First the original pulse splits into two, 
a direct outgoing one and another one going through the 
origin. Both pulses travel along null geodesies taking away 
part of the energy of the excitation toward future null in- 
finity, J^"*", with the help of the massless Higgs field. It 
is found, however, that this way the excited monopole re- 
leases only about half of the energy received. The rest of 
the energy of the original pulse seems to be restrained by 
the monopole in accordance of which it develops a long 
lasting quasi-stable 'breathing state' in the central region 
and certain expanding shell structures in the far away re- 
gion. The frequency characterizing the breathing state 
varies in time and it approaches asymptotically the value 
of the vector boson mass from below. In the far away re- 
gion, where the Yang-Mills and the Higgs fields are prac- 
tically decoupled, the massless Higgs field does the boring 
job of transporting the energy released gradually by the 
central monopole to , while the behavior of the mas- 
sive Yang-Mills field in the far away asymptotic region 
can be characterized by the formation of certain expand- 
ing shell structures where all the shells are built up by very 
high frequency oscillations. These oscillations are found 
to be modulated by the product of a simple time decaying 
factor of the form t~^^'^ and of an essentially self-similar 
expansion. 

The time decay of the examined quantities shows cer- 
tain type of universality. The total energy associated with 



the hyperboloidal hypersurfaces decays in time with power 
—2/3, while the amplitude of the oscillating fields decay 
with power —5/6. In a recent work by Forgacs and Volkov 
Q, based on the use of a linear approximation of the BPS 
monopole, explanation is provided for these universalities. 

The investigated dynamical magnetic monopole is de- 
scribed as a coupled 5(7(2) YMH system. The Yang-Mills 
field is represented by an su(2)-valued vector potential Aa 
and the associated 2-form field Fab reads as 

Fab^\7aAb-\/bAa+ig[Aa,At] (1) 

where [ , ] denotes the product in su(2) and g stands 
for the gauge coupling constant. The Higgs field (in 
the adjoint representation) is given by an su(2)-valued 
function ip while its gauge covariant derivative reads as 
Va-ip ~ Vai/) + iglAajip]. The dynamics of the investi- 
gated YMH system is determined by the action 

S = J [rr(Fe/F"0 + 2Tr(0,^D"V)] (2) 

where e is the 4-dimensional volume element. 

Our considerations were restricted to spherically sym- 
metric configurations yielded by the 'minimal' dynami- 
cal generalization of the static 't Hooft-Polyakov magnetic 
monopole configurations 0,0 (see also Q)- Accordingly, 
the evolution took place on Minkowski spacetime the line 
element of which, in spherical coordinates {t,r,9,(p), is 
ds^ = dt^ - dr^ - [dO'^ + sm^ 6 d(j?), while the Yang- 
Mills and Higgs fields, in the so called abelian gauge, were 
assumed to posses the form 

= -i [■w{r2{de)a - sin e{d(t)) a} + cos0(d(^)J (3) 
i' = Hr„ (4) 

where the generators {tj} (1=1,2,3) of su(2) are related 
to the Pauli matrices aj sts = ^a^ , moreover, w and H 
were assumed to be smooth functions of t and r. 

The field equations relevant for this system are 

r^drW — r^dfw = w [(to^ — l) + g^r^H^^ (5) 
r^d'^H + 2rdrH - r'^d'iH = 2w'^H. (6) 
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The only known analytic solution to (O and ^ is the 
static BPS monopole 
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where C is an arbitrary positive constant. Since this so- 
lution is stable it was used first to check the efficiency 
of our numerical code. Later we considered the complete 
non-linear evolution of a system yielded by strong impulse 
type excitations of this monopole. All the results below 
concerns the evolution of such an excited BPS monopole. 

The only scale parameter of the above described system 
is the vector boson mass m„ = gH^ , where the limit 
value = limr-,cx) H of the Higgs field can in general 
be shown to be time independent Since in the case 
considered here = C / 0, without loss of generality, 
the parameter choice g = = = 1 can be ensured 
to be satisfied by making use of standard rescalings. 

To have a computational grid covering the full physical 
spacetime - ensuring thereby that the outer grid bound- 
ary will not have an effect on the time evolution - the 
technique of conformal compactification, along with the 
hyperboloidal initial value problem, was used. This way 
it was possible to study the asymptotic behavior of the 
fields close to future null infinity, as well as, the inner 
region for considerably long physical time intervals. 

The conformal transformation we used is a slight modi- 
fication of the static hyperboloidal conformal transforma- 
tion applied by Moncrief 01 . It is defined by introducing 
first the new coordinates T and R instead of t and r as 



T = ijt- ^/uj-^r^ + 1 and ^^X^^ll±l — \^ 



(8) 



where u: is an arbitrary positive constant. The Minkowski 
spacetime is covered by the coordinate domain given by 
the inequalities — oo < T < +oo and < i? < 1. Then the 
conformally rescaled metric can be given as Ijab = ^^gab, 
where the conformal factor is f2 — u{l — B?) /2. The 
R = 1 coordinate line represents through which the 
conformally rescaled metric gab smoothly extends to the 
coordinate domain with R> \. 

Using the substitution ^/(f, r) = h{t,r)/r + H^ the field 
equations, lO and ij^J, in the new coordinates read as 

^w = w\^{w''-l)+g^{h + H^Rn-'-fj (9) 
^h^2{h + H^Rn-^)w^, (10) 
where the differential operator is defined as 
4i?2 
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These equations can be put into the form of a first order 
strongly hyperbolic system Q. The initial value problem 
for such a system is known to be well-posed [T(il |. In par- 
ticular, we solved this first order system numerically by 
making use of the 'method of line' in a fourth order Runge- 
Kutta scheme following the recipes proposed by Gustafs- 
son et al 10']. All the details related to the numerical 



approach, including representations of derivatives, treat- 
ment of the grid boundaries are to be published in The 
convergence tests justified that our code provides a fourth 
order representation of the selected evolution equations. 
Moreover, the monitoring of the energy conservation and 
the preservation of the constraint equations, along with 
the coincidence between the field values which can be de- 
duced by making use of the Green's function and by the 
adaptation of our numerical code to the case of massive 
Klein-Gordon fields, made it apparent that the phenom- 
ena described below have to be, in fact, physical properties 
of the magnetic monopoles. 

In each of the numerical simulations initial data on the 
r = hypersurface was specified for the system of our first 
order evolution equations. In particular, a superposition 
of the data associated with the BPS monopole, 0, and 
of an additional pulse of the form 
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with a > b > 0, which is a smooth function of com- 
pact support, was used. This choice, providing non-zero 
time derivative for w, corresponds to "hitting" the static 
monopole configuration between two concentric shells at 
r = a — b and r = a + b with a bell shape distribution. 
Basically the same type of evolution occurs when instead 
of (9tw)o we prescribe {dTh)o in a similar fashion. 

For the sake of brevity, all the simulations shown below 
refer to the same pulse 11211 corresponding to the choice 
of the parameters a = 2, 6 = 1.5, c = 70, d = 10 and 
Lj = 0.05. The energy of this pulse is 55.4% of the static 
monopole. Hence, the yielded dynamical system cannot 
be considered as being merely a simple perturbation of the 
static monopole. We would also like to emphasize that the 
figures shown below are typical in the sense that for a wide 
range of the parameters characterizing the exciting pulse, 
qualitatively, and in certain cases even quantitatively, the 
same type of responses are produced by the monopole 

To start off consider first the product of the energy 
density e and 47rr^ plotted on succeeding T — const hy- 
persurfaces, providing thereby a spacetime picture of the 
time evolution (see Fig.Q. The use of 47rr^e, instead of 
e, makes it easier to see the main characteristics of the 
dynamics up to J^"*". In the early part the direct energy 
transport to by the Higgs field, with the velocity of 
light, is apparent. The developing 'breathing state' of the 
monopole and the formation of the expanding shells of 
high frequency oscillations are both clearly manifested. 

The energy radiated to J^'^ can be pictured by plot- 
ting (see Fig.lJJ the product of the energy current density 
S and 47rr^ against time T at J^^ {R ~ 1)- It is apparent 
that the arrival of the two pulses is followed by a small 
scale but systematic energy loss of the system which has 
exactly the same characteristic period as the inner breath- 
ing state of the monopole. From the logarithmic plot the 
asymptotic behavior 'inr^ Sasympt ~ CgT~'''s can be read 
off where 7^ « 1.66. In virtue of the energy conservation 
this relation implies that the energy associated with the 
T = const hypersurfaces approaches to the energy of the 
asymptotic final state as T~'^^'^, in agreement with Q. 

By inspection of the evolution of the field variables w 
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FIG. 1: Spacetime diagram showing the time evolution of the 
energy density associated with shells of radius r, i.e. of inr'^e. 




FIG. 2: The time dependence of the product of the energy 
current density S and 47rr-^ at ^+ is shown on the intervals 
0.8 < T < 4 and 0.8 < T < 250. 



and h it is obvious that the expanding high frequency os- 
cillations are associated with the massive Yang-Mills field 
exclusively. Fig.^lshows a constant time slice at T = 1.695 
of the evolution of w — wq. The formation of the shells 
built up by high frequency oscillations is transparent on 
Fig.O which is reminiscent of Figs. Land 2. of The 
behavior of these oscillations in the as ymp totic region can 
be explained by referring to results of where this phe- 
nomenon has already been found to characterize the evo- 
lution of the simplest linear massive scalar field. In doing 
this, notice first that by 10 w/r satisfies asymptotically 
the massive Klein-Gordon equation. This, along with the 
claims of Q|, implies then that the oscillations are mod- 
ulated by two factors. By an overall factor t~^^^ scaling 
down the oscillations in time, moreover, by an essentially 
self-similar expansion, i.e. by a function depending on t 
and r only in the combination p — r jt. 

Probably, the most interesting unexpected feature of 
the time evolution is the appearance of the breathing state 
of the monopole (see the central region of Fig.jQ. To have 
a quantitative characterization of this phenomenon it is 



FIG. 3: The difference w — wq is plotted on time slice T = 
1.695. The formation of the high frequency oscillations sorted 
in certain shell structures is transparent. 



informative to consider a constant R slice of the deviation 
e — £o of the full energy density e and that of the static 
monopole eo. The breathing state of the monopole can 
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FIG. 4: The deviation e — eo of the full energy density e from 
that of the static monopole eo is plotted against the physical 
time t for intervals < t < 100 and 100 < t < 6700 at R = 
0.0254. It is striking how quickly the breathing state develops. 

be characterized by the time dependence of the amplitude 
and the frequency of the associated quasi-normal oscilla- 
tions, as well as, by the power spectrum of the oscillations. 
The time dependence of the amplitude is shown by Fig.|3| 
Note that the center of the oscillations is actually lower 
than £o, which implies that the time average of the energy 
contained in the central region is smaller than the energy 
contained in the same region of the static monopole. This 
is most likely due to nonlinear effects. 

The time dependence of the frequency and the ampli- 
tude of the oscillations was determined by fitting a sim- 
ple function of the form e — eo = asm{ujt -I- 6) -I- c to 
the numerical data on a sufficiently short time interval 
so that this interval was shifted point by point through 
the entire time evolution. The resulted graph is shown 
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by Fig.0 The frequency of the oscillations is essentially 
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FIG. 5: The frequency difference niw — u) and the amplitude 
a of the oscillations of the breathing state, shown by Fig.^ is 
plotted against the physical time t. 

increasing, asymptotically approaching the value of the 
vector boson mass rriw The logarithmic plot suggests 
a relatively simple form for the overall behavior of the 
frequency-time function. In particular, the approximate 
relation uj ~ rriw — C^t seems to be valid not merely 
asymptotically but for the entire evolution, where 7^ was 
found up to a high accuracy to be 2/3. Exactly the same 
value 7^ =2/3 was found in Q. The asymptotic behavior 
of the amplitude of the oscillations can be approximated 
by the simple form aasympt ~ Cat , where 7^ ~ 0.833, 
which is in good agreement with the value of —5/6 of @|. 
The energy contained in some finite radius is described by 
the time average of the energy, i.e. by the function c, not 
by the amplitude a. For this reason there is no contra- 
diction between the exponents —5/6 in the central region 
and —2/3 at infinity. The function c appears to be smaller 
and decays faster than the amplitude, and can also take 
positive and negative values depending on the location. 

It is also of interest to consider the power spectrum 
P{uj;ti,t2) of the oscillations defined as twice of the ab- 
solute value of the Fourier transform {e — eo){i^;ti,t2) = 
(271)-'/^//^" (e-eo)(t)e-""'df. Fig.H shows P{Lu-ti,t2) 
with t2 having the fixed value t2 = 11107 and ti being cho- 
sen to take the values 6.8, 42, 60, 95 and 165, respectively. 
By varying ti it is possible to monitor the change of the 
frequency of the oscillations in the relevant early period. 



These graphs support the perturbative result of Forgacs 
and Volkov jF| claiming that there has to be an infinite 
family of 'resonant states', associated with the breath- 
ing state of the monopole, possessing discrete frequencies. 
See, for instance, the peaks of the graph of P{uj; ti = 6.8) 
at w « 0.844, 0.933, 0.9645, 0.9775, 0.985, ... . According to 
Fig.Elthe lower frequency members of the resonant states 
die out faster than the higher frequency ones which is con- 
sistent with the above mentioned increase of the overall 
frequency of the breathing state. 

In summary, the time evolution of strongly excited BPS 
monopoles has been studied. It is found to be generic that 
in the central region a long lasting quasi-periodic breath- 
ing state develops. The behavior of all the examined quan- 
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FIG. 6: The power spectrum of the oscillations of e — so at 
R = 0.0254 relevant for the intervals [ti,i2], with t2 = 11107 
while ti taking the values 6.8, 42, 60, 95 and 165 is shown. 



titles justifies the intuitive expectation that in the inner 
region the system settles down to the original static BPS 
monopole, while the self-similarly expanding oscillations 
disperse asymptotically in the far away region. 
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